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Abstract. The origin of Galactic CRs up the knee energy remains unanswered and provides 
strong motivation for the study of 7 -ray sources at energies above 10 TeV. We discuss recent 
results from ground-based 7 -ray Cherenkov imaging systems at these energies as well as future 
observational efforts in this direction. The exciting results of H.E.S.S. give clues as to the nature 
of Galactic CR accelerators, and suggest that there is a population of Galactic 7 -ray sources 
with emission extending beyond 10 TeV. A dedicated system of Cherenkov imaging telescopes 
optimised for higher energies appears to be a promising way to study the multi-TeV 7 -ray sky. 


1. Introduction 

The origin of hadronic Cosmic-Rays (CRs) remains an unsolved problem since their discovery 
nearly 100 years ago. A major area of focus concerns the Galactic component of CRs, which 
is thought to be dominant at least up to the knee (~ 10^^ eV; 1 PeV), the energy at which a 
spectral break is seen in the CR spectrum. It is generally believed that shell-type supernova 
remnants (SNRs) are responsible for the acceleration of the Galactic component of CRs based 
on energy budget considerations (of SNRs and of the Galactic CRs), and the good agreement 
between the observed properties of the CR spectrum and diffusive shock acceleration (DSA) 
theory. While DSA theory is able to explain many aspects of the observed Galactic CRs, a 
controversial aspect concerns the maximum particle energies attained ED H. A maximum 
energy limit of ~ 10^^ eV is generally predicted m which falls somewhat short of the knee 
energy. Considerably higher energies approaching the ankle (~ 10^^ eV) are thought possible 
if one accounts for signihcant fluctuations in the magnetic held [H El Alternative source 
scenarios have also been considered such as local Gamma-Ray-Bursts (CRB) |,'I11117j . In the 
conventional context of SNR acceleration of CRs, so-called Superbubbles, which beneht from the 
combined effects of many SNR shocks and possibly also from the stellar winds of OB associations, 
may also accelerate CRs to well beyond the knee dill El- 

We are therefore left with critical questions concerning the origin of CRs at and above the 
knee. Unfortunately, identifying localised sources of Galactic CRs at such energies is impossible 
by direct detection since the CR trajectories are are randomised by Galactic magnetic helds. 
However, the interaction of GRs with local matter can produce a hux of 7 -rays^ in the multi-GeV 


^ and neutrinos 


to multi-TeV energy band which act as an accessible tracer of CR accelerators. Moreover, the 
energy spectrum of the 7 -ray flux is expected to follow that of the parent CR spectrum at the 
site of 7 -ray production, thus providing important information concerning the CR acceleration 
properties. We note that the interaction of CRs with ambient matter proceeds with a typical 
inelasticity of order ~0.15. A source of 7 -rays, when interpreted in the above framework, 
is therefore a source of hadronic CRs of energies roughly a factor 10 higher. Establishing 
beyond doubt a hadronic origin for a 7 -ray flux is non-trivial since we must often deal with 
additional leptonic sources of 7 -ray production such as inverse-Compton scattering and non- 
thermal Bremsstrahlung. These latter components are often diminished at higher energies due 
to the strong cooling suffered by electrons in magnetised post-shock environments. Morphology 
differences between the hadronic and leptonic components are also expected, due to a number 
of physics issues such as energy-dependent diffusion. Spatially-resolved spectra and energy- 
resolved morphology studies of 7 -ray sources are necessary to disentangle the hadronic and 
leptonic components. In particular, establishing the spectra and morphology of sources of 7 -rays 
at energies 10 TeV and above is vital to understanding where and how hadronic CR acceleration 
to the knee is taking place. 

The detection of TeV 7 -rays is the exclusive domain of ground-based methods, and a major 
goal of the field is to establish the types of sources responsible for Galactic CRs. Given the 
maximum energy problem discussed above, it is not sufficient however just to establish numerous 
7 -ray source types and their spectral slope, but to clearly establish the maximum energies of their 
7 -ray fluxes well into the multi-TeV domain. The performance of ground-based 7 -ray detectors 
at energies E > 10 TeV will therefore attract increasing interest over the coming years. 

2. Present Ground-Based Results at Energies E > 10 TeV 

The stereoscopic imaging atmospheric Cherenkov technique is employed by all of the present 
ground-based detectors such as H.E.S.S., CANGAROO-III, VERITAS and MAGIC (stereoscopy 
will be employed with an additional telescope - MAGIC-II). These systems have built on 
the experience from earlier detectors like the HEGRA lACT-System UHl which established 
stereoscopy, and the Whipple 10 metre telescope m, which first demonstrated the power 
of imaging. Survey instruments such as the particle detector arrays MILAGRO, TIBET, 
KASGADE and CASA-MIA for example also operate in the TeV to PeV domain. We will 
however not cover results from these experiments, a selection of which may be found in 
uni nni na m cni and references therein. 

The TeV results from H.E.S.S. (see summary in m) have signalled a turning point for 
TeV 7 -ray astronomy in that the number of sources has grown rapidly, and individual source 
morphology is revealed for the first time [SIEl- For the shell-type SNR RX J1713.7—3946 there 
are now sufficient photon statistics to permit energy-resolved morphology, and spatially-resolved 
spectral studies HO]. Over half of the more than 20 TeV 7 -ray source have a Galactic origin^ [Sj. 
Many of these new Galactic sources, with integral fluxes in the range few % to 15% Grab (at 
E > 1 TeV), have plausible associations with nearby shell-type SNR or plerionic SNR (pulsar- 
wind-nebulae). A striking observation is that the majority of the new sources exhibit generally 
hard photon spectra, typically characterised by an unbroken power-law dN/dE E ^ with 
index T in the range 2.0 to 2.5, and extending to up to 40 TeV. The majority are also extended 
in morphology, with diameters ranging from a few arcmin to a degree. These results are in 
line with expectations of DSA theory and give strong hints as to the type of Galactic GR 
accelerators. The establishment of spectra above energies ~10 TeV is often limited however by 
statistics, given the available observation time and sensitivity of H.E.S.S. Deeper observations 
will in fact be performed on selected H.E.S.S. sources such as RX J1713.7—3946, which will be 

^ This is somewhat biased by the H.E.S.S. observation programme, which is split roughly 50:50 between Galactic 
and extragalactic targets. 


observed at high zenith angles (this method is discussed shortly) during 2005 in a strong effort 
to establish the spectrum beyond 40 TeV. 

H.E.S.S. is optimised for 7 -ray detection in the energy range ~ 100 GeV up to a few 10s of 
TeV. The specification of the four H.E.S.S. telescopes, with mirror area 120 m^, 0 . 12 ° camera 
pixel resolution, 5 deg camera field of view (EoV), and a telescope spacing of 120 m, gives 
H.E.S.S. an optimal sensitivity in the 1 to 10 TeV range, where the CR background rejection 
is strongest, and the angular resolution is better than a few arcminutes (see sensitivity curve 
Figure ^). In this energy range an energy flux sensitivity of ~10“^^ erg cm“^s“^, or 1% Crab 
flux, in 25 h observation time at low zenith angles (LZA) is achieved. The performance at 
higher energies begins to decrease since the camera FoV, and telescope spacing (defining the 
physical size of the array) imposes a limit on the detection efficiency. This limitation may be 
offset somewhat by employing high-zenith angle (HZA) observations [30] , where the observation 
zenith angle z is larger than ~ 55°. At such angles, the longer atmospheric path length means 
that the air shower is viewed at greater distances from the telescope. As a result, higher energy 
EAS preferentially trigger the detector from considerably larger distances (impact parameters) 
than they would at lower zenith angles, thereby increasing the effective collection area for a 
given energy, but at the same time increasing the energy threshold. Using this technique, the 
effective collection area of a ground-based system like H.E.S.S. at high energies approaches 1 km^ 
(Figure 121). 

HZA observations have already been demonstrated effectively in the study of the Crab Nebula 
by the HEGRA lACT-System [ 2 ] and Mkn 421 by H.E.S.S. [21 . The energy spectrum of the 
Crab is now established up to energies approaching ~70 TeV after very deep observations 
(~380 h) with the HEGRA lACT-System (Fig. jS)). The mechanism put forward to explain 
the (unpulsed) TeV emission of the Crab centres on the inverse-Compton upscattering of soft 
photon fields by accelerated electrons. The inverse-Compton 7 -rays seen up to 70 TeV, along 
with the synchrotron emission seen up to GeV energies by EGRET strongly imply the presence 
of electron acceleration to PeV energies | 2 ]. Thus, the Crab is likely the first source established 
as a cosmic Pevatron. It therefore remains to be seen if it, and other Galactic sources are 
capable of accelerating hadrons to similar energies. H.E.S.S. observations of Mkn 421 were 
triggered by a flaring state reaching fluxes of up to ~6 Crab. At the H.E.S.S. site this northern 
source is visible at zenith angles ~ 65°, yielding an energy threshold about 2 TeV, however 
with collection areas exceeding 1 km^. Despite the fact that the energy spectrum of Mkn 421 
exhibits an exponential cutoff {dN/dE E ^exp(E/Ec) where E^ = 2> TeV) limiting the 
photon statistics above 10 TeV, HZA observations permitted a detection rate of up to 8 7 -rays 
per minute. For the first time, variability studies on time scales of hours were performed on this 
source in the multi-TeV energy range. The HZA observation method has limitations however. 
Observability issues come into play as most sources do not spend much time at high zenith 
angles. Indeed this method is better suited, for instruments in the south, to viewing sources 
with northern declinations in order to maximise the available observation time per night (the 
opposite of course is true for southern sources viewed from the north). Other issues such as 
atmospheric stability can also reduce the overall efficiency of the observations. 

3. A Dedicated Ground-Based Instrument for Energies E > 10 TeV 

It is clear that we are only just beginning to sample the 7 -ray sky at energies ~ 10 to 100 TeV 
at 1% Crab flux sensitivties. Over the past decade, is has become apparent that the stereoscopic 
atmospheric Cherenkov imaging technique has proved highly successful for 7 -ray detection in the 
100 GeV to ~10 TeV energy range. It is therefore worthwile investigating the same technique 
optimised for higher energies. The Monte-Carlo simulations of Plyasheshnikov et al. m 
first looked at the potential of the stereoscopic atmospheric Chereknov imaging technique for 
E > 10 TeV 7 -ray astronomy, and we draw on those results here. They considered a cell of four 
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Figure 1. Energy flux sensitivities for various instruments operating in the 7 -ray domain, 
compared the measured and predicted (synchrotron and inverse-Compton) fluxes from the Crab 
Nebula. For the ground-based instruments, limited here to Cherenkov imaging systems “100- 
GeV Arrays” (representing H.E.S.S., VERITAS and CANGAROO-III), HEGRA lACT-System, 
Whipple, 5@5), a 50 h 5 it detection limit is required. See ini for details of the MAGIC sensitivity. 
For the space-based detectors (GLAST, INTEGRAL, MEGA) a 10® sec obs time is specified. 
The reported 7 -ray fluxes from HEGRA, GELESTE, EGRET and COMPTEL are also shown. 
From 1^. See Figure 01 for a more recent Crab spectrum from HEGRA. 



Figure 2. H.E.S.S. effective collect¬ 
ing area vs. energy (following stan¬ 
dard 7 -ray selection cuts) for a vari¬ 
ety of zenith angles (20°, 45°, 65°) 
jS]. Vertical lines indicate the en¬ 
ergy threshold (peak differential rate 
for a Grab-like spectrum) applicable 
to each zenith angle. An increase in 
the effective collecting area at higher 
energies is clearly noticed for higher 
zenith angle observations. 
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Figure 3. Energy spectrum 

of the Crab Nebula from HEGRA 
lACT-System observations [ 2 ]. Data 
from low and high zenith angles are 
compared. The spectrum extends to 
~70 TeV. The highest energy spectral 
point has an excess significance +2.7 cj. 


imaging Cherenkov telescopes, equipped with focal plane cameras of photomultiplier’s (PMTs) 
of various sizes and telescope spacing. The main motivation was to investigate the potential of 
the technique, when applying generally conservative specifications to the telescope design. In 
this sense a quite conservative quantum efficiency (QE) (photons to photoelectrons) value of 0.1 
was assumed. Monte-Carlo simulations of both 7 -rays and protons viewed at the vertical were 
used. Two different pixel sizes for the focal plane detectors (or cameras) were considered, 0.3° 
and 0.5°, and for each case, a two-level tailcut scheme to remove pixels dominated by skynoise 
was also applied^. A telescope trigger was achieved when at least 2 two adjacent pixels meet 
a defined threshold (10 to 20 p.e.). A stereoscopic trigger was demanded when at least two 
telescopes have triggered from the same EAS. In order to gauge the performance characteristics 
of the cell, a moment-based image parametrisation scheme according to Hillas was applied. This 
employed the width parameter and its stereoscopic counterpart mean-scaled-width IT P as the 
image shape basis for CR background rejection, while the image major-axes were used in the 
reconstruction of arrival directions. A point spread function (PSE) of size ~0.15° to 0.1° (std. 
deviation of a fitted Gaussian) was indicated for the cell when assuming a pixel spacing of 0.3°. 
These trigger conditions and analysis algorithms applied to the cell were based on the methods 
used in the HEGRA lACT-System, and now, also employed by H.E.S.S. 

It is worth highlighting here the major conclusions of the simulation study: 

• Telescope Mirror Area: Eor energies E > 10 TeV, only a modest mirror area A ~ 10 m^ 
is required to detect sufficient Cherenkov photons (or equivalent photoelectrons). Eigure0 
depicts for both 7 -ray and CR proton induced air-showers, the number of Cherenkov 
photoelectrons (p.e.) contained within an image vs. impact parameter (the distance of 
the air-shower core to the telescope) and energy. Experience has shown that ~50 to 100 
p.e. per image are required to accurately reconstruct individual image parameters such as 
width, length and orientation. A mirror area of 10 m^ is therefore more than sufficient to 

® This is also referred to as image cleaning. In H.E.S.S. data analysis, tailcuts of 5 and 10 p.e. are often used. 













Figure 4. Average photoelectron 
yield for Cherenkov images vs. impact 
parameter (R = distance to the air 
shower core) for a range of energies. 
Solid curves ( 7 -rays); dashed curves 
(CR protons). A field of view of 
infinite diameter is used, to accept 
photons from all directions. From 

m- 


detect and reconstruct events with energy F > 10 TeV out to impact parameters 500 m or 
more. 

• Focal Plane Camera FoV: EAS that trigger a telescope out to 500 m will be imaged with 
considerable angular displacement in the focal plane, especially when observed at low zenith 
angles (0° to 45°). Figure [3 shows the relationship of the image centroid distance from the 
focal plane centre (the so-called dist parameter) vs. impact parameter. Images with impact 
parameters approaching 500 m will trigger up to 4° off-axis. When one allows for edge 
effects of the finite FoV, a typical image length of order 0.8° and the fact that many sources 
might be extended, fields of view in excess of 8 ° diameter would be considered necessary. 
Depending on the pixel sizes (0.3° or 0.5°), one would require between 300 to 800 pixels in 
order to achieve such fields of view. 

• Telescope Spacing: The efficiency (the fraction of events triggering relative to those 
incident) of the stereoscopic trigger depends on the distance between telescopes, on the 
telescope FoV, and energy of the EAS. The results of Figure 0 suggest that a large telescope 
spacing of order 500 m is suited to multi-TeV energies. Figure El depicts the detection 
efficiency for four telescopes with spacing 500 m. In fact a spacing in the range 300 to 
500 m was found necessary to achieve an 80% efficiency over the desired energy range, 
when employed in combination with very large FoV cameras. 

Assuming then a telescope spacing of 500 m, an energy threshold (the peak differential rate) 
of 7 to 10 TeV was achieved for a point source with a Crab-like energy spectrum. For a source 
with a Crab-like flux, detection rates of 3.5 7 -rays and ~700 CR background events per hour 
(over the entire FoV for the CR background) were achieved. A combination of a cut on the 
angular distance 0 < 0 . 1 ° between the reconstructed and assumed directions in addition to a 
cut on the scaled-width parameter (eg. W < 1.0) accepted <0.03% of CR background events, 
whilst retaining 50% of 7 -ray events. After these cuts, a Crab-like source yields detection rates 
of 1.7 7 -rays per hour against a CR background rate of 0.2 events per hour, giving a signal 
to noise ratio approaching a factor 10 and excess significance ~ 2cT/\/hr^. It is interesting to 

^ Here we assume that the CR background is estimated from a region in the FoV with solid angle a factor 10 
larger than that of the source. Eq. 17 from Li & Ma 1281 with a = 0.1 is then used to estimate the significance. 





Figure 5. Average value of dist (the image 
centroid distance from the focal plane centre) 
vs. impact parameter R. Images are viewed 
and reconstructed using a large FoV camera 
(9.3°, with pixel sizes 0.3°) and pixel trigger 
threshold of 15 p.e. Prom m 
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Figure 6. The efficiency (fraction of events 
triggering relative to those incident) of event 
trigger vs. the detector FoV for several 7 -ray 
energies. A telescope spacing of 500 m is used 
and pixel trigger threshold of 10 p.e. From 
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compare these numbers with the detection rate obtained on real Crab data, for energies F > 7 to 
10 TeV from both the HEGRA lACT-System and H.E.S.S. detectors. For HEGRA, a detection 
rate of ~2 7 -rays and 0.3 GR per hour is achieved at 2(T/\/hr [2]. For H.E.S.S., the Crab 
culminates at middle-range zenith angles of ~ 45°, improving the collecting area somewhat at 
higher energies (FigureEJ. On the Crab, H.E.S.S. achieves a 7 -ray rate near culmination of ~ 6 
7 -rays and 0.4 CRs per hour at ~ 5(T/\/hT- If one factors in the improved collection area for the 
H.E.S.S. Crab result, the performance of a single cell of small telescopes is quite similar to that 
of H.E.S.S. and the HEGRA lACT-System at the high energies considered. Another important 
performance measure is the energy resolution, which is estimated at ~25%. This is sufficient to 
reconstruct energy spectra accurately and in particular determine cutoff features. 

Of course the majority of sources under future study will emit at sub-Grab flux levels. The 
H.E.S.S. results suggesting a population of sources in the 5 to 15% Crab flux range (F > 1 TeV), 
would demand considerable observation time from a single cell. This is helped somewhat by the 
hard spectra of the new sources (which would have fluxes ~ 0.2 Crab-flux level for F > 10 TeV). 
In order to detect sources with high statistics in reasonable time (< 100 h), [2Zj also suggest the 
use of multiple cells of 4-telescope arrays. For example they estimate a total of 25 telescopes 
(arranged as 16 cells of 4-telescopes) would achieve detection rates (after cuts) on the Crab of 
~ 1250 7 -rays and 150 CRs in 50 h observation. This observation time would be sufficient for 
extraction of a spectrum for a 0.1 Crab-flux source. 

It should again be stressed that the performance of a high energy array of telescopes would be 
improved for hard-spectra sources (the simulation study assumed a softer Crab-like spectrum), 
like those being uncovered by H.E.S.S., and indeed those expected to be the sources of hadronic 
CRs. Finally we note that the detectability of 7 -ray sources will be affected by the attenuation 
of 7 -rays due to pair production on the microwave background m- This loss strongly peaks at 





~1 PeV with a mean free path of ~10 kpc, meaning that PeV 7 -ray astronomy would be limited 
to searching for relatively nearby Galactic sources. 


4. Conclusions 

Studies of 7 -ray sources at energies above 10 TeV are vital if we are to understand how and 
where Galactic cosmic-rays are accelerated to and above the knee. Since the present ground- 
based telescopes such as H.E.S.S. are optimised for energies E <10 TeV, it is suggested that a new 
dedicated system of telescopes be employed to open up the 7 -ray sky at higher energies. Monte- 
Gar lo simulations of a cell of four modest-size telescopes employing the same stereoscopic imaging 
technique as in H.E.S.S. and the HEGRA lAGT-System suggest that sufficient sensitivity is 
achieved in the E > 10 TeV range. These encouraging results, from rather conservative 
simulations, suggest that the stereoscopic imaging approach is a viable method to explore the 
multi-TeV sky. More detailed simulations are now underway. 
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